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Abstract 

Given a positive energy solution of the Klein-Gordon equation, 
the motion of the free, spinless, relativistic particle is described in a 
fixed Lorentz frame by a Markov diffusion process with non-constant 
diffusion coefficient. Proper time is an increasing stochastic process 
and we derive a probabilistic generalization of the equation (dr)^ = 
— ^dXiidX^. A random time-change transformation provides the 
bridge between the t and the r domain. In the r domain, we obtain 
an M'^-valued Markov process with singular and constant diffusion 
coefficient. The square modulus of the Klein- Gordon solution is an 
invariant, non integrable density for this Markov process. It satisfies 
a relativistically covariant continuity equation. 

PACS number: 03.65.Bz 

Running title: Free relativistic particle 
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1 Introduction 



The purpose of this paper is to discuss a possible scenario for the free, spin- 
less, relativistic particle, extending on Nelson's stochastic mechanics. This 
challenging problem has generated considerable interest in the past twenty 
years, see [|], ^, |], ||, 0, |14|, |T8|, [T^, ^ ^ , and references therein. As is well- 
known, the main difficulty in this generalization stems from the non existence 
of nontrivial Markov diffusion processes possessing the required relativistic 
covariance, [0, |12|. A number of attempts have been made to circumvent this 
problem, for instance by considering different classes of stochastic processes, 
see e.g. |,|,[Tg,0. 

In this paper we propose a new approach. Corresponding to a positive 
energy solution of the Klein-Gordon equation, motion of the particle is de- 
scribed in a fixed Lorentz frame by an M^-valued Markov diffusion process 
with non-constant diffusion coefficient (local covariance matrix). The role 
of proper time is played here by an increasing stochastic process, namely 
a quadratic variation process. We also derive a natural probabilistic gen- 
eralization of the equation (dr)"^ = —^dX^dX^. A random time-change 
transformation provides the bridge between the t and the r domain. In the r 
domain, we obtain an M^-valued Markov process with singular and constant 
diffusion coefficient. The fourth component is proportional to a stopping 
time. The square modulus of the Klein-Gordon solution is an invariant, 
non integrable density for this Markov process. It satisfies a relativistically 
covariant continuity equation. 
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2 Background on diffusion processes 



In this section, we review some essential concepts and results of the kine- 



matics of stochastic mechanics. We refer to [22|- |24], M for a thorough 



account. In order to avoid any confusion, stochastic processes will be de- 
noted by capital letters, as it is customary in probability. Let {fl,£,P) be a 
probability space. A stochastic process {X(t);to < t < ti} mapping [to, ^i] 
into L'^i^, S, P) is called a finite- energy diffusion with constant diffusion co- 
efficient /nCr^ if the path X{uj) belongs a.s. to C([to, ti]; M") (n- dimensional 
continuous functions) and 



X{t)-X{s)= J' p{T)dT + a[W+{t)-W+{s)], to<s<t<ti, (2.1) 

where the forward drift f3{t) is at each time t a measurable function of the past 
{X(r);0 < T < t}, and W^{-) is a standard, n-dimensional Wiener process 
with the property that W+(t) — W+{s) is independent of {X(r); < r < s}. 
Moreover, j3 must satisfy the finite-energy condition 

E y^'' p{t) ■ f3{t)dt^ < oo. (2.2) 

In P], Follmer has shown that a finite-energy diffusion also admits a reverse- 
time differential. Namely, there exists a measurable function 7(t) of the 
future {X(r);t < r < ti} called backward drift, and another Wiener process 
W_ such that 

X{t)-X{s)= f\{T)dT + (T[W^{t)-W^{s)], to<s<t<ti. (2.3) 
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Moreover, 7 satisfies 

E lit) ■ i{t)dt'j < 00, (2.4) 

and W_(t) — W_{s) is independent of {X(r); t < r < ti} . Let us agree tfiat 
dt always indicates a strictly positive variable. For any function / defined 
on [to, ti], let 

d^fit) ■.= f(t + dt)-f{t) 
be the forward increment at time t, and 

d.f{t) = f{t)-f{t-dt) 

be the backward increment at time t. For a finite-energy diffusion, Follmer 
has also shown in that the forward and backward drifts may be obtained 
as Nelson's conditional derivatives, namely 

and 

the limits being taken in L^(f2,i3, P). It was finally shown in |^ that the 
one-time probability density p{-,t) of X{t) (which exists for every t > to) is 
absolutely continuous on M" and the following relation holds Vt > 

E{P{t)-i{t)\Xit)} = a^Vlogp(X(t),t). (2.7) 

Let us introduce the current drift v{t) := (/3(t) + 7(t))/2 and the osmotic 
drift u{t) := (Pit) — 7(t))/2. Notice that, when a tends to zero, v tends to 



X, and u tends to zero. The finite-energy diffusion X(-) is called Markovian 
if there exist two measurable functions &+(■,■) and &_(•,■) such that = 
b+{X(t),t) a.s. and ■y(t) = b^{X(t),t) a.s., for all t in [to,^i]- The duality 
relation ( p.7| ) now reads 

6+ (X (t) , t) - 6_ (X (t) ,t) = a^V log p(X(t) , t) . (2.8) 

This immediately gives the osmotic equation 

2 

«(x,t) = yVlogp(x,t), (2.9) 

where M(a;,t) := {b^{x,t) — b_{x,t))/2. The probability density p(-, ■) of X(t) 
satisfies (at least weakly) the Fokker-Planck equation 

dp „ . ^ 



The latter can also be rewritten, in view of (|2.8|) , as the equation oj continuity 
of hydrodynamics 

^ + V-(t;p) = 0, (2.10) 

where v{x, t) := (6+(x, t) + 6_(x, t))/2. 

Nelson's stochastic mechanics ||2^ , |^ is a quantization procedure for 
classical dynamical systems based on diffusion processes. Given a quantum 
evolution {■?/'(x, t\t<^ < t < t\\ ^ namely a solution of the Schrodinger equation 

satisfying Carlen's finite action condition 

l|V^||^GLL[[to,+oo)], (2.12) 



it is possible to construct a measure on path-space under which the canonical 
coordinate process X is a finite-energy Markov diffusion process. Writing 
ipix, t) = exp R{x, t) + ^S{x, t), we have that the process X has current and 
osmotic drift fields given, respectively, by v{x,t) = ■^VS{x,t), and u{x,t) = 
—'VR{x,t). In particular, the (forward) Ito differential of X is given by 



dX{t) 



-V (3fJ log^(X(t), t) + !3 log ij{X{t),t)) 
m 



dt 



'-dWJt), (2.13) 
m 



see [0], |T|, Chapter IV], and references therein. Moreover, the probability 
density p{-,t) of x(t) satisfies 

p{x,t) = \ij{x,t)\', VtG[to,ti]. (2.14) 

We now need to recall the generalization to Markov processes with diffusion 
coefficient of the form /„o"(x,t)^, with a{x,t) G R, the above kinematics, see 
e.g. |2T[, [gg, ^ and jg. Let {X{t);to < t < ti} be an ra-dimensional 



Markov diffusion process whose increments admit the two representations 

d+X{t) =b+{X{t),t)dt + a+{X{t),t)d+W+{t), (2.15) 

d_X{t) =b4X{t),t)dt + a_{X{t),t)d^W4t), (2.16) 

where the drifts and diffusion coefficients are sufficiently regular. Here, 
and W- are standard, n-dimensional Wiener processes adapted to the past 
and the future, respectively, of the process X. Then, the following relations 
hold 

(j+(x, t) = (T_(x, t) = (t(x, t) (2.17) 

b+(x, t) - 6_(x, t) = -V (a(x, tfp(x, t)) . (2.18) 
P ^ ' 
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Letting, as usual, v := (6+ + &-)/2 denote the current drift, we have, as 
before, that the Fokker-Planck equation can be rewritten as the continuity 
equation (|2.1CI| ). Let R{x,t) := ^ log (cr(x, t)^p(x, t)). Then, observing that 



(|2.18| ) can be rewritten as 

b+{x,t) - b^{x,t) = a{x,tyV\og (a{x,ty p{x,t) 
we get the following expression for the osmotic drift u 

u{x, t) = cx{x, t fVR{x, t). (2.19) 

3 Description in a fixed inertial frame 

For the definitions and results about continuous martingales that occur in this 
section, we refer the reader to the Appendix and, for a thorough treatment, 
to [0 and |T^. Let ip := exp[i? + ^S] be a solution of the Klein-Gordon 
equation 

on [0, oo), such that S satisfies on the same time interval 

— + cVm^c^ + V5 ■ V5 = 0. (3.2) 

ot 



This is precisely the family of solutions that was considered in [|T9| . Indeed, 
let 

() = ^ {f^V - f^v) (x, t) = \f{x, t)\' X (^VSfe ()) . (3.4) 



First of all, observe that because of (|3.1|), the pair (p, j) satisfies a continuity 
type equation 

Moreover, from (|3.2|) it also follows that 

P > 0, (3.5) 
j-j-cV = -c2M^ (3.6) 



The latter two conditions are equivalent to (5.10) in [[191 , p. 4705]. This class of 
solutions is non empty. It contains at least all positive-frequency plane-wave 
solutions, namely solutions of the form 

ip{x,t) = exp[{i/h){p ■ X - p'^t)], 

where 



p° = +cJ p ■ p + m'^c'^. 



Let 



It follows from ( |3.2| ) that po{x) > 0. Let us suppose that po integrates to one. 
Suppose that {X{t);t > 0} is a Markov diffusion process such that X(0) is 
distributed according to po, and having forward differential 

(3.7) 
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where is a standard, 3-dimensional Wiener process. In the previous 

notation, we have 

(J^(x,t) := _j_gs( ,\- 

A standard calculation shows that p{x,t) in ( |3.3D satisfies the Fokker- 
Planck equation corresponding to ( p.7|) . Hence, p{x,t) is the probability 
density of X{t). Moreover, the current drift of X is 

/ ±\/s \ 

t;(X(t),t):=pV^ (X(t),t), (3.8) 



mc^ dt 



and the osmotic drift is given by 



uiXit),t)=(^^)iXit),t). (3.9) 

Let us define 

M\t):= I d^Wlis), t>0 z = 1,2,3. (3.10) 

The stochastic processes {M*(t); t > 0} are continuous /oca/ martingales with 
the same quadratic variation 

< M > (*) =< M' > it) ■= [ _'^ ds. (3.11) 

mc^ as ^ \ / ' / 

Moreover, their cross variations satisfy 



< M\ > (t) = 0, 1 < z 7^ J < 3, Vt > 0. (3.12) 

Observe that < M > (t) is a strictly increasing process with differentiable 
sample paths. It satisfies, in view of ( p.2|) and (|3.8|) , 

d<M > it) = -—^^ — -—-rdt 
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\ m?c^ + VS{X{t),t) ■ VS{X{t),t) 



dt 



_ VS{X{t),t)-VS{X{t),t) 
\\ " " w?c^ + V5(X(t), t) ■ VS{X{t),t) 



1 - 



dt 



v{X{t),t)-v{X{t),t) 



dt. 



(3.13) 



If for some t and u; we have v{X{t,uj),t) = 0, then d < M > (t) = dt. We 
can therefore think of d < M > {t) as of an increment of the particle random 
proper time. Notice that ( |3.13| ) imphes 

d<M > (t) 



dt 



< 1, a.s. 



Hence, 



< M > (t) <t, a.s., Wt > 0. 



(3.14) 



Furthermore, we have the following property oi d < M > (t). Let us first 
recall that, in view of (OI) and (12.61), we have 



(3.15) 



where the symmetric increment of X at time t is defined by 

d+X{t) + d^X{t) X{t + dt) -X{t-dt) 
d,X{t) :- . 

We can therefore rewrite ( |3.13| ) as follows. 



^[d<M> {t)X 



[icf + v{X{t),t) ■ v{X{t),t) {dt) 



,d{ict) 2 



dt 



Y + v{X{t),t)-v{X{t),t) 



{dtf 



{d{ict)f + E{dsX{t)\X{t)} ■ E{d,X{t)\X{t)} + o{{dtf 



(3.16) 
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Let us introduce 



X{t) 
ict 



Then ( p.l6| ) can be written as follows 



[d<M> {t)f = --E{d,X,it)\X,it)} ■ E{dMt)\X,{t)} + oiidtf), 

(3.17) 

generalizing the relation (rfr)^ = —^dX^dX,, of classical relativistic mechan- 



ics. 



4 From the t to the r formulation: A random 
time- change 

An important consequence of Levy's characterization of the Wiener process, 



see e. g. |16, p. 82], is that any continuous local martingale may be viewed 
as a time-changed Wiener process, cf. Theorem A. 2 in the Appendix. For 
r > 0, let us introduce the stopping time 

T(r) := inf{(T > :< M > (a) > r}. 

(In the case when the probability of the event {uj : lim(_^oo < M > {t,uj) = 
oo} is strictly less than one, the stopping time T(r) has to be suitably mod- 
ified, see [|T3|, pp.174-175]). Notice that 

<M> (T(r)) = r, r(< M > (t)) = t. (4.1) 

Then, the processes 

W^l(r) :=M^(,), ^ = 1,2,3, 
12 



are, standard, one- dimensional Wiener processes. Moreover, in view of prop- 
erty (p.l2|) , we can apply a theorem of F. B. Knight (see Theorem A. 3 in 
the Appendix), and conclude that the processes VF|(r) are pairwise inde- 
pendent, and can, consequently, be viewed as the components of a standard, 
three-dimensional Wiener process W+{t). Let us introduce the stochastic 
process 

X(r) :=X(T(r)), r > 0. 
In view of (|3.7|) , (|3.11|) , and ( ^!T|) , the forward differential of X is given by 



d+X{T) = -VS + -VR (X(r), T{T))dT + J-d+W+ir). (4.2) 



m 



m 



m 



Notice that X{t) is a non-Mar kovian M'^- valued diffusion process with con- 
stant diffusion coefficient ^Is- The Markov property has therefore been 
destroyed by the random time change. From the first relation in ([4.1|) , we 

get 



c/T(r) 
dr 

In view of ( p. Ill ), we then get 

rfr(r) 



d < M > 
Jt 



1 -1 



(T(r)) 



' ^^-(X(r),T(r)). 



dr mc^ dt 

Let us now define X^{t) := icT^r), S{x^) = 5'(x,a;^) := S{x,t), and R{x^) 
R{x,x^) := R{x,t). Equations (|4.2| ) and ( [4.3|) can now be rewritten as 



(4.3) 



d+X{T) 
d^XUr) 



^-V~S + -VR] {X{r),X\r))dr + 



1 dS 



m 



m 



dx^ 



{X{T),X\T))dr. 



-d+W^{T) (4.4) 



(4.5) 
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Notice that Xu{t) defined by 



r 



X(r) 
XJt) 



is an M^- valued Markovian stochastic process with marginal density at r = 
given by po{x)-6{x^) since T(0) = a.s. Let us introduce the diffusion matrix 



1 



\ 







Then ( [4 .41) - ([4 .51) can be rewritten in the form 



m 



+ S^V,i?j (X,(r))dr + Sd+W^,(r), (4.6) 



where Wi, is any standard, four- dimensional Wiener process whose first three 
components form W+. The Fokker-Plank equation is then 



dfhy 
dr 



(4.7) 



Because of (PTTD, ^p{x^) := exp[R{x,y) + j-S{x„)] satisfies 



-if. 



The latter is equivalent to the system of p.d.e.'s 



2m 2m 



mc 



Let 



m 2m 



0. 



(4.8) 
(4.9) 
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From ([4.9|), we get that p(xy) satisfies ( |4.7| ). Thus, iJi{dxy) := p{xi,)dxi, is an 



invariant (a-finite) measure for (|4.4|) - (|4.5| ). Notice that p{x,^) also satisfies 
the manifestly covariant equation 

V. ■ (p-Vj) = 0. 

m 

If we give p{dx^) := p{xy)dxi, as initial measure to (^^), the probabilistic 
picture is lost. Nevertheless, it is possible to make sense of the time-reversed 
diffusion along the lines of pp. 44-45]. Moreover, consider the generator oi 
Xy acting on smooth functions with compact support 

m J 1 



The adjoint of this operator with respect to the measure p{dxy) p. 104] 
is given by 

\m / 2 

Hence, in spite of the lack of a probabilistic picture we can still think of 
;^Vy5 as of a bilateral velocity field associated to the equilibrium measure 
p{xy). Then, recalling that ( |3.2| ) implies 

5 



Vi,S{xi,) ■ VyS{xy) = -rn^c 



we get a r-domain counterpart of relation (p.l7|). 



5 Conclusion and outlook 



Although the process ( [4.41) - ( ^TSf ) is not relativistically covariant, our descrip- 
tion does not appear to be in confiict with classical relativistic mechanics nor 
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with Nelson's nonrelativistic theory. Indeed, in the case when we let the noise 
intensity in ( p.7|) tend to zero, we recover the equations of classical special 
relativity. In particular, the quadratic variation < M > defined by ( p.llj ) 
tends to proper time r and T(r) tends to ^(t). On the other hand, when 
V ■ V « c^, motion of the particle in a fixed inertial frame is described by 
a Markov diffusion process with diffusion coefficient nearly equal to ^. The 
Schrodinger equation is then recovered through some appropriate procedure 
neglecting terms multiplied by /i^/c, cf. |1T9| , p. 4708]. 

Acknowledgments: The author wishes to thank Paolo Dai Pra for some 
helpful comments on the results of this paper, and Lorenza Viola for several 
useful conversations on relativistic stochastic mechanics. 
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A Continuous Martingales 

We collect in this appendix a few basic facts about continuous martingales. 
We refer the reader to [jl5| and [^ for the proofs and more information. 
Let {fl,A,P) be a probability space, and let := (J-'t), t G [tc^i], be a 
non- decreasing family of sub a -algebras of A. It will be always assumed 
that J-'tQ contains all the zero-probability sets in A, and that the filtration 
is right-continuous, namely {nJ^t,t > s} = Ts-, for every s. A stochastic 
process X := {X{t);t G [toj^i]} is said to be continuous if it has continuous 
trajectories X{t, u),t G [to, ti], with probability one. It is called J-'- adapted ii, 
for every t in [toj^i]? ^{t) is J^rineasurable. In that case, it is customary to 
write X = {X(t), J^t)t(^[to,ti]- The stochastic process M = {M{t), J^t)te[to,ti] is 
called a martingale (with respect to the filtration J-'t,t G [to7^i]) if it satisfies 
the two following conditions: 

1. E{\M{t)\} < oo, Vt, 

2. E{M{t)\J^s} = M{s), t > s. 

The continuous martingale X = {X(t),J^t)telto,ti] is called square-integrable 
if E{X(ty} < oo, Wt in [to7^i]- Let X = {X(t), J^t)te[to,ti] be a continuous, 
square-integrable martingale with X(to) = with probability one. Then, 
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the celebrated Doob-Meyer decomposition theorem imphes that there exists 
a unique representation for X{t)^ 

X{tf=<x>{t) + M{t), te[to,ti], 

where < X > is an adapted, continuous, increasing process with < X > 
(to) = 0, and M is a continuous martingale with M(to) = 0. The process 
< X > {■) is called the quadratic variation of X{-). Consider now two 
continuous, square-integrable martingales X — {X{t), J^t)te[to,ti] ^^'^ ^ ~ 
{Y{t), J^t)telto,ti]- Their cross-variation process < X,Y > is defined for t e 
[to:ti\ by 

< X,Y > (t) :^^[< X + Y > (t)- < X - Y > (t)] . 

The cross variation of X and Y is characterized by the fact that XY— < 
X, y > is a continuous martingale. 

A nonnegative random variable T{u!) is called a Markov time or a stopping 
time (relative to the filtration J^t,t e [to, ti]) if, for all t in [to, ti], 

{cu : T{uj) <t}e Tf 

Let T be a stopping time relative to Tt,t e [io,^i]- Then, the cr-field Tt 
of events determined prior to the stopping time T consists of those events 
A e ^ for which [A D {a; : T{ui) < t}] G J^t, for every t in [to, ti]. 

The stochastic process M = {M(t), J-'t)t>o, is called a local martingale 
if there exists an increasing sequence of stopping times (with respect to JF) 
{Tn}, n — 1,2, . . . such that 

1. P(limT„ = oo) = 1; 
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2. for any n, the "stopped process" {M(t A Tn),t > 0, is a martingale. 

The decomposition result for the product of two square-integrable martingale 
extends to local martigales as follows. Let X = {X(t),J^t)tGlto,ti] and Y = 
(Y(t), J^t)te[to,ti] be two continuous, local martingales. Then there is a unique 
adapted process < X,Y > such that XY— < X, F > is a continuous local 
martingale. We write < X > instead of < X,X >. We can now state P. 
Levy's martingale characterization of the Wiener process. 

Theorem A.l Let M = {M{t), J-'t)t>o be a continuous local martingale with 
M(0) = a.s. and quadratic variation < M > (t) = t. Then M is a standard 
Wiener process. 

As a corollary to this theorem, we have the following result showing that con- 
tinuous local martingales can be viewed as time-changed Wiener processes. 

Theorem A. 2 Let M = {M{t), J-'t)t>o be a continuous local martingale with 
M(0) =0 a.s. Suppose that limf^oo < M > (t) = oo with probability one. 
Define, for r > 0, the stopping time T{t) := infjt > :< M > (t) > r} . 
Then the time-changed process W = {W{t) := M(T(r)),^^ := J^t{t))t>o is 
a standard Wiener process and we have, with probability one, 

M{t) = W{< M > (t)), t > 0. 

The multivariate extension of this theorem is due to F. B. Knight [^, p. 179]. 
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Theorem A.3 Let M = {M{t) = {M^{t), . . . , M"(t)), :Fj)t>o} be a contin- 
uous process with M(0) = a.s. Suppose that the components M* are local 
martingales satisfying the two following conditions: 

1. limt_^oo < > if) — i — 1, . . . ,n a.s.; 

2. < M\ > (i) = 0, i 7^ J, t > 0. 

Let T,{t) := m{{t > :< M' > (t) > r},t > 0,i = 1, . . . ,n. Then the 
W^ir) := M\Ti{T)),T > 0,i = l,...,n are independent, standard Wiener 
processes. 

The nontrivial content of this theorem is that, ahhough the M* are not 
independent, applying first the appropriate time- changes, and then forgetting 
the time changes, we get independent Wiener processes. Forgetting the time- 
changes consists in replacing the filtrations {Q^} with the poorer filtrations 
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